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We study cascaded harmonic generation of hybrid surface plasmons in integrated planar waveg¬ 
uides composed of a graphene layer and a doped-semiconductor slab. We derive a comprehensive 
model of cascaded third harmonic generation through phase-matched nonlinear interaction of fun¬ 
damental, second harmonic and third harmonic plasmonic modes supported by the structure. We 
show that hybrid graphene-semiconductor waveguides can simultaneously phase-match these three 
interacting harmonics, increasing the total third-harmonic output by a factor of 5 compared to the 
non-cascaded regime. 
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I. INTRODUCTION 

Parametric wavelength conversion in general and har¬ 
monic generation in particular are among the most in¬ 
tensively studied phenomena in optics and photonics^. 
It was demonstrated that important advantages can be 
obtained by cascading several parametric processes^’^. 
For instance, since a wide range of materials provide 
stronger quadratic nonlinear response compared to cubic 
nonlinear response, the efficiency of the third harmonic 
generation (THG) can be strongly enhanced by utilis¬ 
ing a cascaded process consisting of second harmonic 
generation (SHG) and sum-frequency generation (SFG), 
which mixes fundamental wave (FW) and second har¬ 
monic (SH) to generate the third harmonic (TH). The 
conversion efficiency in this regime is the highest when 
the phase-matching (PM) conditions for both SHG and 
SFG are satisfied. 

Doped graphene is an atomically thin tunable plas¬ 
monic material^’^, which can be incorporated into vari¬ 
ous components in nanoscale and microscale optics. In 
particular, graphene demonstrates strong optical nonlin¬ 
earity, especially when utilised near the plasmon reso¬ 
nance®^®. It is therefore of significant interest to establish 
a theoretical framework and explore the viability of cas¬ 
caded nonlinear optical processes in graphene-enhanced 
wavelength conversion devices. 

In this work we develop a theoretical model for cas¬ 
caded third harmonic generation in graphene-based plas¬ 
monic waveguides. We show that engineering modal dis¬ 
persion allows to simultaneously satisfy SHG and SFG 
phase-matching for efficient cascaded THG, and offers 
5 times increase of TH output compared to direct non- 
cascaded THG. We consider a planar waveguiding system 
consisting of a graphene sheet placed over a doped semi¬ 
conductor sandwiched between two dielectric layers, see 
Fig. 1. The dielectric permittivity of a doped semicon¬ 
ductor is described by the Drude model with the plasma 
frequency comparable to the Fermi energy in graphene. 
This geometry allows the coexistence of three plasmonic 
modes®, which as we show below can be phase-matched, 
strongly increasing effective nonlinear optical interaction. 


In our model we assume that the nonlinear interaction of 
propagating plasmons is achieved due to the nonlinear 
response in graphene. In this work we focus only on 
nonlinear processes in the waveguide; the way of linear 
in-coupling and out-coupling of light to and from such 
hybrid-graphene waveguides was studied elsewhere^®. 

This paper is structured as follows. In Sec. H, we 
derive the expressions for the nonlinear conductivities 
at FW, SH and TH in graphene. In Sec. HI, we de¬ 
rive the SHG and THG coupled-mode equations and 
the corresponding nonlinear optical coefhcients, includ¬ 
ing direct and cascaded contributions to the TH. In 
Sec. IV, we study the phase-matching and present the 
results of a numerical simulation, showing that hybrid 
graphene-semiconductor waveguides can simultaneously 
phase-match these three interacting harmonics, increas¬ 
ing the total third-harmonic output by a factor of 5 com¬ 
pared to non-cascaded regime. 
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Figure 1. (Color online) Geometry of a planar hybrid 
graphene-semiconductor waveguide. A waveguide consists of 
a monolayer graphene sheet placed at the height d 2 above the 
doped semiconductor with permittivity Sp and thickness di 
sandwiched between two dielectric layers with permittivities 
£i and £ 2 . The third dielectric layer with permittivity £3 is 
placed on top of a graphene layer. The optical pump at the 
frequency uj (red arrow) is generating the second harmonic at 
the frequency 2 a; (green arrow) and the third harmonic at the 
frequency 3a; (blue arrow). 
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II. NONLINEAR CONDUCTIVITIES OF 
GRAPHENE 


First, to characterize the nonlinear response of a 
graphene layer, we methodologically derive expressions 
for nonlinear conductivities taking into account in-plane 
plasmon momentum, which results in nonzero quadratic 
conductivity and enables SHG. To describe nonlinear 
properties of doped graphene, we employ a quasi-classical 
description based on the solution of the Boltzmann equa- 
tion^^“^^. It is applicable at low frequencies, fno < Ep, 
where Sp is the Fermi energy. In the collisionless limit, 
the electron distribution function /(r,p,t) satisfies the 
kinetic equation written as 


df , df 


dr 


+ ’^p ® ® ^ ) b:; “ ® ’ (1) 


df 


dp 


dE 

where for Dirac’s massless quasi-particles Vp = — 

op 

the energy £ = vpp, vp ^ c/300 is the Fermi 


velocity, and e = —|e| is the charge. Following the per¬ 
turbation approach^’we solve this equation using iter¬ 
ations. 


Here we consider the case of p-plasmons or TM- 
polarized incident wave with the tangential electric field 
component of the form = Exoe*'^“ (and qvpjuj <C I). 
In this section, we assume that a graphene monolayer is 
placed in the plane xy. We note that there is no compo¬ 
nent of magnetic held normal in respect to the graphene 
surface. Because of the graphene nonlinearity, the held 
excites higher harmonics. We look for the distribu¬ 
tion function in the form 


/ = foie) + + {fid + + c.c. 


( 2 ) 


hrst, the second, and the third order oscillating correc¬ 
tions (the superscripts are omitted below), respectively, 
whereas the tangential electric held is as follows 

E^ = -k c.c. + -k c.c. 

(3) 

+ c.c. + ... . 


The electric current density at w is given by the hrst 
order correction to the Fermi-Dirac distribution function 

= 4e^ Vp^(p) = 4ey y , 

where 

/„= (4) 

i(UJ — Vp^q) luj \ oj J p dp 

Therefore, the hrst order current amplitude can be writ¬ 
ten as follows 
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4e^g dfo 
iuj dp 


vpcos'^Q (5) 




This is a known result^^. In the local semiclassical limit, 
the linear conductivity of graphene can be reduced to the 
Drude model, which only takes into account intraband 
processes 


^intra(^) 


ie^ £p 
pffl LO ’ 


It can also be corrected to include relaxation 


ie^ £p 
~ Fh? (w + *Ti-\ J 


( 6 ) 

(7) 


where /o(£) is the Fermi-Dirac distribution (step-like at 
low temperatures, kpT «C £f), fiH^ /sw 


where is the relaxation rate. 

The second order correction is dehned as follows 
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It gives rise to the electric current density at 2uj j 2 u> = turns zero. Integrating the second term leads to the fol- 
4eX]''^p/2aj(p) ■ Integrating the hrst term in Eq. (8) re- lowing expression 
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which can be reduced to 


3 e^vl 






8 nhPuj^ 

Therefore we obtain the nonlinear source 


( 10 ) 


3 e^vl 


J2uj — —Xq- 


8 Trh?uj^ 


qE'^ exp(—i2a;t + i2qx) , (11) 


which is in agreement with the earlier results^^s^ Thus, 
considering frequency conversion with plasmons, the 
quadratic conductivity of graphene cr 2 {uj) can be defined 
as follows 


0 - 2 ( 0 ;) 
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Next, the third order correction is defined as 

^ df^\ 

J3ui — XT 7 -r I eEi^— - \-eE2ujF .— I 
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To obtain the contribution to the third harmonic current 
= 4eX) Vp/ 3 tj(p) from the first term 
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where q = ksp{ijj) is the plasmon wavenumber supported 

by a waveguide at ui. we expand this term as follows 
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Here integrating the expressions proportional to q^ re¬ 
turns zero. After integrating the terms proportional to 
q^ and the current amplitude can be written as follows 

J 3 LI = = [^ 3 ( 0 ;, (? = 0) -b Acr3(a;, g)] , (16) 

where the third order conductivity, in the homogeneous 
case^^“^^d6 given by 


sponding to the second term given by 


f(ii) _ 
J 3uj 


1 p % 

3i{(jj - Vp^q)^ dpx 


1 


eE< 


2ui 


JL 

dpa 


{-i} 

\%UJ \ 




-3i:;o + 

p up 


('+¥) 


(19) 


and get the x-projection of the nonlinear current as fol¬ 
lows 


( 73 ( 0 ;, g = 0) 


.1 / evF\^ £f .1 Vp 

8 irh? J cu 8 nh'^ Spuj^ ’ 

(17) 


is corrected by an additional term 


= \^2{u2)E2^E^ , (20) 

which entirely corresponds to the cascaded process in 
THG. 


Acr3(a;,g) 


. 5 vf / qvF\‘ 
16 wh? £pLo^ V w / 
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(18) 


Further, by analogy with the second harmonic contri¬ 
bution f 2 uj, we derive the third harmonic current corre- 


III. COUPLED MODE EQUATIONS 

Neglecting losses, the dispersion relation and 
p-polarized eigenmodes of a hybrid graphene- 
semiconductor waveguide were analysed in Ref.®. 
Similar to the approach outlined in Ref.^®, to describe 
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the plasmon conversion due to the graphene nonlinear¬ 
ity, we derive equations for the modal envelopes. The 
frequency conversion can be modeled by coupled equa¬ 
tions for the amplitudes of FW, SH and TH plasmonic 
modes, Ai, A2 and ^3, respectively. Assuming that the 
conversion efficiency is small, we utilize the following 
magnetic field ansatz 

Hy{x,z,t) = 

+ e-*2-M2(/rz)/i2(2w, (21) 

where /ii.2,3{a;) are the transverse profiles of the modes. 
Disregarding the influence of the generated modes on the 
pump, the coupled-mode equations can be written as fol¬ 
lows 

—h 7 iAi = 0 , (22a) 

^ + 72A2 = , (22b) 

^ + 73^3 = g2AiA2e-*['=?p(3-)-fe=p(2-)-fe=V(-)l^ (22c) 
dz 

_|_ ^3y[3g-d*:ap(3“)-3fe)p(pj)]z ^ 

where the term proportional to 52 is responsible for the 
so-called cascaded contribution to the THG^’^. Assuming 
that the FW plasmon is launched at z = 0 , we study 
the plasmon frequency conversion modeled by Eqs.( 22 ). 
The conversion from the TH to the SH can be neglected, 
because, as will be shown below, the TH efficiency in 
this kind of graphene-semiconductor waveguides remains 
much smaller than the SH efficiency, Asiz) <C A2(^). 

Loss coefficients 71,2,3 in Eqs.( 22 ) can be obtained 
perturbatively from the dispersion equation, introduc¬ 
ing dissipative corrections, namely the linear frequency- 
dependent surface conductivity of graphene 

a{uj) = a^^^Lo) + ( 23 ) 


and the dielectric permittivity of the semiconductor 


ep{uj) = 1 - 


U}„ 


uj(uj + iv, 


_ 1 


IVU}„ 




( 24 ) 


where is a damping coefficient. The coefficients g2, §2, 
(73 in nonlinear sources can be derived from the equations 
for the energy flows in the corresponding modes as follows 




Figure 2. (Color online) (a) Dispersion of three supported 
plasmon modes in a hybrid graphene-semiconductor waveg¬ 
uide (solid blue lines), corresponding asymptotes (dashed 
black lines), corresponding points of phase-matching for FW, 
SH and TH (blue circles); (b) FW, SH and TH mode pro¬ 
files shown for the continuous in-plane electric field compo¬ 
nent. Bold black line corresponds to the graphene layer, hne 
black lines correspond to the interfaces between semiconduc¬ 
tor/dielectric layers. 


Here e^’^(a;) and hy’^{x) are the transverse profiles of the 
SH and TH modes normalized in such a way that at the 
graphene layer ez{d2) = 1 . The field components at the 
p-polarized surface wave are connected as 


ez{.x) 


1 dhy 1 dcx 
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( 25 a) 

( 25 b) 


Finally, we can write the nonlinear coefficients as follows 


92 = - 


^2 

4G2’ 




( 28 ) 


These coefficients define the strength of the nonlinear 
interactions between FW, SH and TH plasmons. Phase 
matching, another important parameter that defines the 
conversion between these harmonics, is discussed in the 
following section. 
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Figure 3. (Color online) Spatial evolution along the 
propagation direction of the normalized SH intensity 
\^ 2 {z)\‘^/\Ai{0)\^ (solid green line), total TH intensity includ¬ 
ing direct THG and cascaded THG \A 3 {z)\^/\Ai{Q)\^ (solid 
blue line), and the direct TH intensity |^3(2)|^/|^i(0)p in 
the absence of the cascading effect assuming ^ 2=0 (dashed 
blue line). 


IV. PHASE-MATCHING AND WAVELENGTH 
CONVERSION 

Next we show that three guided modes corresponding 
to FW, SH and TH can be phase-matched, so that at 
some w, dependent on the geometrical parameters (sepa¬ 
rations between the layers, dielectric permittivity distri¬ 
bution and doping levels), they have almost equal effec¬ 
tive indexes /3^(a;) « /3^(2a;) « f3^{3uj) (or, equivalently, 
klp{uj) = fcgp(2a;)/2 = fcgp(3w)/3). To calculate the dis¬ 
persion and the mode profiles we use the transfer matrix 
method^^. 

To demonstrate this behaviour we choose the follow¬ 
ing set of realistic waveguide parameters. The permit¬ 
tivities of the dielectric layers shown in Fig. 1 are chosen 
as £i = 8 , £2 = 2, and £3 = 1 . For simplicity we do 
not take into account the dispersion in the dielectric lay¬ 
ers. Adding the dielectric dispersion will require adjust¬ 
ing the thicknesses of the dielectric layers to compensate 
for it. We choose the thickness of a semiconductor layer 
as di = 25 nm, and the thickness of the dielectric layer 
between the semiconductor and graphene as ^2 = 20 nm. 
Then for the Fermi energy £p = 0.5 eV, plasma fre¬ 
quency ujp = 0.95fF) relaxation parameter Tintra = 0.3 
ps and the damping coefficient v = O.OOlWp, we obtain 
simultaneous phase-matching for FW, SH and TH, as 
demonstrated in Fig. 2(a), for the fundamental wave¬ 


length in free space of 11.7 ^m. In Fig. 2(b) we show 
the corresponding overlap between the FW, SH and the 
TH modes. Strong graphene nolinearity allows efficient 
nonlinear interaction even between the modes with mod¬ 
erate field confinement, which can be further improved 
by optimizing the waveguide parameters. Remarkably, 
for different dielectric and semiconductor parameters, the 
TH mode can be phase-matched in the backward direc¬ 
tion, although in this paper we focus on the case when all 
three modes are phase-matched in the forward direction. 

In Fig. 3 we demonstrate the evolution of the inter¬ 
acting FW, SH and TH plasmons. We choose a realistic 
pump intensity of 1 MW/cm^ on a graphene layer^®“^^. 
The resulting SHG and THG conversion efficiencies and 
optimal conversion lengths are comparable to the previ¬ 
ous works®d 0 > 22 ^ Particularly, choosing y dimension of 1 
cm and the pump Poynting flux per unit length 1 W/cm, 
the outcome energy flows reach 2 x 10“^ and 3 x 10“® 
W/cm, respectively, at the distance of about 2 fim. We 
note that the direct contribution to THG governed by 
the term proportional to 33 in Eq.(22c) is an order of 
magnitude smaller than the cascaded one. To show this 
explicitly, as a dashed line in Fig. 3, we plot the TH inten¬ 
sity evolution calculated by formally setting g 2 = 0. The 
difference is remarkable, because the plasmon wavevec- 
tor can notably exceed that in free space, leading to the 
efficent SHG process stimulating a cascaded conversion 
to the third harmonic plasmon. 

V. CONCLUSIONS 

To conclude, in this work we have established the the¬ 
oretical framework for the direct and the cascaded THG 
in graphene-semiconductor waveguides, including a de¬ 
tailed derivation of all coefficients in the coupled-mode 
equations based on the waveguide parameters. We have 
shown that these waveguides can support simultaneous 
phase-matching of FW, SH and TH plasmon modes, and 
that the cascaded THG can be several times stronger 
than the direct THG. 
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